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Abstract 

We use ab initio global structural prediction, and 
specifically the minima hopping method, to ex¬ 
plore the periodic table in search of novel oxide 
phases. In total, we study 183 different compo¬ 
sitions of the form MXO 2 , where M=(Cu, Ag, 
Au) and X is an element of the periodic table. 
This set includes the well-known Cu delafossite 
compounds that are, up to now, the best p-type 
transparent conductive oxides known to mankind. 
Our calculations discover 81 stable compositions, 
out of which only 36 are included in available 
databases. Some of these new phases are poten¬ 
tially good candidates for transparent electrodes. 
These results demonstrate, on one hand, how in¬ 
complete is still our knowledge of the phase-space 
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of stable ternary materials. On the other hand, 
we show that structural prediction combined with 
high-throughput approaches is a powerful tool to 
extend that knowledge, paving the way for the ex¬ 
perimental discovery of new materials on a large 
scale. 

Introduction 

Most of the knowledge painfully accumulated over 
centuries concerning the crystal structure of (in¬ 
organic) materials is nowadays gathered in gener¬ 
ally available databases. The most used of these, 
the Inorganic Crystal Structure Database (ICSD) 1 
contains around 170,000 entries. These entries in¬ 
clude unfortunately duplicated structures, insuffi¬ 
ciently characterized phases (e.g., missing the po¬ 
sitions of H atoms), and several alloys. If we re¬ 
strict ourselves to well-defined crystal structures 
we find around 40,000 entries. 

It is this amount of 40,000 entries that has been 
studied theoretically, in a systematic manner, over 
the last few years using high throughput tech¬ 
niques. The results of these large-scale studies can 
be found in excellent publicly available databases, 
such as the Materials Project, 2 the Open Quantum 
Materials Database, 3 or the Ab-initio Electronic 
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Structure Library AFLOWLIB. 4 These are then 
used as a starting point for sophisticated machine 
learning techniques that try to design new materi¬ 
als with tailored properties. These techniques have 
been touted as the most cost effective path to the 
discovery of new materials in a diverse range of 
applications, such as Li batteries,' 5 6 thermoelec¬ 
tricity,® 9 photovoltaics. 10-13 

But let us go a step back, and look at the 40,000 
materials we know nowadays. The first question 
that comes to our mind is if this set is represen¬ 
tative of the number of (thermodynamically sta¬ 
ble) materials that we can create in a lab. The 
answer to this question is probably no. It is true 
that elementary substances and binary compounds 
are relatively well studied, but enormous gaps still 
exist in our knowledge of ternary and multinary 
materials. Of course, a systematic experimental 
endeavor of synthesis and characterization of all 
possible phases is extremely expensive and time- 
consuming. It is for this task that numerical sim¬ 
ulations appear as the most cost-effective way to 
explore the gigantic search space of compositions 
at our disposal. 

In this context, most theoretical studies follow a 
very simple recipe. It is well-known that Nature 
often chooses the same solution to similar prob¬ 
lems. Therefore, one takes a known, experimen¬ 
tally characterized structure, change its chemical 
composition, and hope that this will also be a (dy¬ 
namically and thermodynamically) stable phase 
for the new composition. Sometimes the chemical 
composition can be varied by “brute force” or one 
can use sophisticated machine learning algorithms 
to predict what are the most promising substitu¬ 
tions. 14-17 A recent example of the first approach 
is the work of Carrete and coworkers, 9 who com¬ 
puted all possible compositions of half-Heusler 
compounds in the cubic structure (around 80,000 
possibilities), while the latter approach was used, 
for example, by Hautier and coworkers^ 18 who in¬ 
vestigated Nature’s missing oxides, i.e., which ox¬ 
ides were thermodynamically stable, but still un¬ 
known to mankind. It is evident that in either case 
the major drawback of this approach is the im¬ 
possibility to discover any material with different 
crystal structure than those of the compounds al¬ 
ready contained in available databases. 

In this Article, we address this problem and go a 


step further, showing that structural prediction al¬ 
gorithms can be used on a large scale together with 
high-throughput methods to predict the lowest- 
energy crystal structures of materials with unre¬ 
ported chemical compositions. The problem that 
global structural prediction methods try to solve is 
very simple to explain: given the chemical compo¬ 
sition of a solid, obtain the minimum energy crys¬ 
tal structure, i.e. its ground-state. In possession 
of this structure, we can then evaluate its energy 
and a number of spectroscopic properties using 
standard methods based on density functional the¬ 
ory and beyond. Of course, the number of possi¬ 
ble minima in a solid increases exponentially with 
the size of the unit cell, which makes structural 
prediction an extremely hard numerical problem. 
Fortunately, several smart algorithms for global 
structural prediction appeared over the last decade, 
ranging from random search to evolutionary algo¬ 
rithms.™ 

Our method of choice is the minima hopping 
method. 23 24 This is an efficient algorithm de¬ 
signed to predict the low-energy crystal structures 
of a system given solely its chemical composition. 
At a given pressure, the enthalpy surface is ex¬ 
plored by performing consecutive short molecular 
dynamics escape steps followed by local geome¬ 
try relaxations, taking into account both atomic 
and cell variables. The initial velocities for the 
molecular dynamics trajectories are chosen ap¬ 
proximately along soft mode directions, allowing 
efficient escapes from local minima and aiming 
towards lower energy valleys. Revisiting already 
known structures is avoided by a feedback mech¬ 
anism. The minima hopping method has been 
used for structural prediction in a wide range of 
materials, 25-28 including the dependence on pres- 
surJ 29 and the exploration of binary phase dia¬ 
gram 30 with remarkable results. 

We use this machinery to study a subset of 
oxides, namely those having a composition of 
the type (Cu,Ag,Au)XC> 2 , which includes the “fa¬ 
mous” delafossites Cu(Al,In,Ga)0 2 Pl2ll These 
compounds of this family are still the most promis¬ 
ing p-type transparent conductive oxides (TCO) 
known to mankind. TCOs possess the uncom¬ 
mon property of being at the same time transpar¬ 
ent to the visible spectrum and good electric con¬ 
ductors. These properties make them indispens- 
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able for many high-technology devices which re¬ 
quire transparent contacts, such as flat panel dis¬ 
plays, touch screens, thin-film and stacked so¬ 
lar cells, functional windows, etc. Good electron 
(n-)doped TCOs, namely those based on SnC> 2 , 
In 2 C> 3 , and ZnO, are already widely used in com¬ 
mercial applications. However, potential p-type 
TCOs identified up to date have conductivities 
at least one-two orders of magnitude lower than 
their n-type counterparts and carrier mobilities too 
small for large-scale exploitation. The best p- 
type TCO until now is CuCr ( j_ x .,Mg v 02, 35 which 
displays a conductivity of 220 0 _1 cm _1 and a 
hole mobility of about 1 cm 2 /Vs, while also suf¬ 
fering from poor transparency, with transmission 
in the visible smaller than 30%. The origin of 
the relatively higher hole mobility of Cu delafos- 
sites relies on the fact that the highest valence 
bands are obtained through the strong hybridiza¬ 
tion of almost-degenerate oxygen 2 p and copper 
3d states. This hybridization reduces the local¬ 
ization of the top valence states on oxygen atoms, 
leading to more dispersive p-d anti-bonding bands 
with smaller hole effective masses. 36 A tetrahe¬ 
dral coordination of oxide ions (as in delafossite 
crystals) is particularly advantageous as it allows 
strong hybridization. Cu 1+ (or equivalently Ag 1+ 
and Au 1+ ) appear ideal elements for creating a p- 
d dispersive top valence while preserving trans¬ 
parency, as their closed d shell will prevent from 
absorption in the visible. Even if it is believed that 
a Cu +1 configuration is the best to obtain TCOs, 
we decided nevertheless to explore blindly also 
compositions that favor Cu 2+ or even Cu 3+ , as 
we do not want to preclude the possible formation 
of crystalline structures different from delafossite 
that would let emerge interesting electronic prop¬ 
erties in different environments. 

We can extract from the analysis above some 
expected good rules for the design of improved 
p-type TCOs: (i) cations should have d shells 
proximate in energy to oxygen 2 p states; (ii) 
cation d shells should be closed to avoid opti¬ 
cal absorption in the visible range; (iii) strong 
hybridization of oxygen p and cation d states is 
required to increase the band dispersion. 36 Fol¬ 
lowing these simple ideas, few other Cu oxides 
were already successfully tested in experiments af¬ 
ter Cu(Al,In,Ga)02J 31 ' 34 such as CuCrMg02, 3 ^ 


SrCu202^ or LaCuOS^ and (Cu,Ag)Sc0 2 . 1 ^ 
Our aim is now to extend and make more sys¬ 
tematic this investigation by pre-screening in sil- 
ico all possible Cu, Ag, and Au based ternary ox¬ 
ides, in order to offer to experimentalists a reli¬ 
able guide on the stability and electronic properties 
of the still unknown compositions. While we se¬ 
lected the broad family of (Cu,Ag,Au)X02 oxides 
as some of their members are examples of com¬ 
pounds satisfying the empirical rules for good p- 
type TCO candidates, we decided at the same time 
to limit as much as possible the restrictions im¬ 
posed to the choice of chemical elements inside 
this family, not to hinder the emergence of unre¬ 
ported crystal structures with completely different 
electronic properties. 

Note that the set under consideration here is 
also a subset of the systems studied using high- 
throughput techniques in the above-mentioned 
work of Hautier and coworkers. 18 As a conse¬ 
quence, we will be able to provide a direct com¬ 
parison of the success rate of our mixed ap¬ 
proach, that combines high-throughput with struc¬ 
tural prediction, in comparison with a purely high- 
throughput study. For reference, we should keep 
in mind that three metastable materials of the 
type (Cu,Ag,Au)X02 unreported in experimental 
databases were obtained in Ref. [T8l namely (i) a 
trigonal phase of AgCo02 (50 meV above the con¬ 
vex hull of thermodynamic stability) based on the 
prototype experimental structure of Agln 02 , (ii) a 
tetragonal phase of AgCs02 (340 meV above the 
hull) based on the prototype CuCs 02 , and (iii) an 
hexagonal phase of AgLa 02 based on the proto¬ 
type AgA102 (4 meV above the hull). 

Methodology 

We performed minima hopping simulations for 
all stoichiometries of the type (Cu,Ag,Au)X02, 
where X is any element of the periodic table up 
to Bi with the exclusion of the rare gases and the 
lanthanides. Forces and energies were calculated 
within density functional theory 4041 in the projec¬ 
tor augmented wave (PAW) formalism 42 as imple¬ 
mented in VASP. 43 44 For a given stoichiometry, 
the initial geometries were obtained randomly, en¬ 
suring only that the minimal distance between the 
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atoms was at least equal to the sum of the covalent 
radii. The MHM searches were performed using 
the Perdew-Burke-Erzernhof (PBE) 45 approxima¬ 
tion to the exchange-correlation functional. We 
used default “high” accuracy energy cutoffs. Each 
minima hopping run was repeated at least twice, 
using both one and two formula units (4 or 8 
atoms). This may seem a relatively small number 
of atoms in the unit cell, but if we look at the exper¬ 
imentally known structures of the (Cu,Ag,Au)XC >2 
family, we realize that almost all of them have a 
ground-state with less than 8 atoms (exceptions are 
some metastable phases of AgCC >2 with 16 atoms 
per unit cell and a stable phase of AgBC >2 with 128 
atoms per unit cell). Therefore, we assume that 
this is not a major limitation. 

We then compared the structures obtained in our 
runs with the ones present in available experimen¬ 
tal and theoretical databases. 23 Almost all exper¬ 
imental structures appeared at the early stages of 
our minima hopping runs, which certainly proves 
the efficiency of the method for this kind of task. 
Finally, we took the experimental structures and 
other relevant theoretical phases that we discov¬ 
ered during our minima hopping runs, and used 
them as prototypes for modified stoichiometries. 
We believe that this procedure, which effectively 
mixes high-throughput techniques with structural 
prediction is essential in such large scale appli¬ 
cations as each method can be used to provide a 
check on the other. 

In total we investigated 183 stoichiometries, and 
we obtained 1,000 minima, of which ~7,500 
were further analyzed. In this last step we fol¬ 
lowed the same protocol as in the Materials Project 
database: spin-polarized calculation using the 
PBE 45 exchange-correlation functional, with the 
exception of the oxides of Co, Cr, Fe, Mn, Mo, 
Ni, V, W where an on-site Coulomb repulsive in¬ 
teraction U® with a value of 3.32, 3.7, 5.3, 3.9, 
4,38,6.2, 3.25, and 6.2 eV, respectively, was added 
to correct the d- states. PAW setups were taken 
from the version 5.2 of VASP. At this stage the en¬ 
ergy cutoff was set to 520 eV (irrespective of the 
elements considered) and k-point grids were au¬ 
tomatically chosen to ensure convergence to bet¬ 
ter than 2 meV per atom. For all low-lying min¬ 
ima we studied the thermodynamic phase equi¬ 
libria of the ternary system, considering the en¬ 


ergy balance with respect to all possible decom¬ 
positions in ternary, binary and elementary com¬ 
pounds that respect the overall stoichiometry: i.e, 
we measured the thermodynamic stability by cal¬ 
culating the energy distance from the convex hull 
of stability, which is the set of lines that con¬ 
nects the lowest energy ordered phases, using PY- 
matgen 47 and the data available in the Materi¬ 
als Project database. According to this definition, 
a compound is stable if its total energy distance 
to the convex hull is zero. Finally, the crystallo¬ 
graphic analysis of the structures was performed 
using FINDSYM.® 

Results 

Our results are summarized in Fig. [l] In these pe¬ 
riodic tables the color scale indicates the distance 
of the calculated energy of the ground-state struc¬ 
ture to the convex hull of stability. Light green 
cases indicate stable structures (i.e. that the spe¬ 
cific composition correspond to a crystal structure 
on the hull). We further indicate the atomic num¬ 
ber (upper left), the space group we predict for the 
lowest energy phase (upper right), and the space 
groups of the lowest-energy experimental structure 
that we found in available databases (in parenthe¬ 
sis below the chemical symbol). This means that 
all squares without a number below the chemical 
symbol represent crystal structures predicted theo¬ 
retically. Further details can be found in the Sup¬ 
porting Information that includes structural and 
electronic data of the most important (meta-) sta¬ 
ble structures. 

As we can see from Fig. [TJ many compositions 
are thermodynamically stable, and many others are 
quite close to the convex hull. Note that our cal¬ 
culations are performed at zero temperature and 
pressure for perfect periodic crystals. However it 
is well known that unstable phases can be stabi¬ 
lized by temperature, pressure, defects, dopants, 
etc. Moreover, due to the theoretical error associ¬ 
ated to the calculation of the total energy, an in¬ 
version of the ordering of the phases very close 
in energy is always possible. In view of the above, 
we believe it is relevant to discuss all compositions 
that are either thermodynamically stable or quasi¬ 
stable (within 20 meV from the convex hull). This 
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Figure 1: Distance to the convex hull (in meV per atom) for Q 1 XO 2 , AgXC> 2 , and A 11 XO 2 . The colors 
indicate the distance to the convex hull of stability, light green meaning that the composition is thermody¬ 
namically stable. For each element we show the symbol (center), atomic number (top left), space group of 
the lowest energy structure (top right), and the space groups of the phases found in experimental databases 
(bottom, in parenthesis). 
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choice of 20 meV as a threshold of stability comes 
from the observation that well studied experimen¬ 
tal compounds such as CuInC^ and AgNiC >2 are 
above the hull by 10-20 meV in our calculations 
performed for perfect stoichiometric bulk crystals. 
The set of structures that we select within this 
threshold of 20 meV contains therefore phases that 
have a large chance of being synthesized experi¬ 
mentally. We find that this stability condition is 
fulfilled by 81 compositions (30 containing Cu, 29 
containing Ag, and 25 containing Au), of which 
only 36 (19 containing Cu, 11 containing Ag, 
and 8 containing Au) are present in experimental 
databases. We are now in the position to compare 
our predictions with previous theoretical works. 
We compared directly the crystal structures when¬ 
ever they were available for comparison. If the 
structure in literature was not given, we used avail¬ 
able information on the symmetries and distance 
to the convex hull to compare with our results. 
There are two phases (AgLaC >2 and AgCs02) pre¬ 
viously predicted in Ref. 18, while theoretical 
crystal structures for CuBa02, CuAuC> 2 , CuHgC> 2 , 
AgLi02, AgNaC >2 can be found in the Materials 
Project database. In the two cases of AgLaC >2 
and AgNaC >2 we find the same crystal structure 
as previously predicted. However, in the other 
five cases the MHM runs were able to find lower- 
energy phases. In particular, the reported struc¬ 
tures of CuBa02, AgLi02, and AgCsC >2 were well 
above the hull, while our lowest-energy structures 
are thermodynamically stable. We should also ob¬ 
serve that there may exist even lower energy struc¬ 
tures that were not detected in the MHM runs (be¬ 
cause their unit cell contains more than 8 atoms, 
for example). However, the compounds we predict 
to be thermodynamically stable have a high chance 
of being experimentally synthesizable, even if the 
actual atomic arrangement may eventually differ 
from the one that we have determined. In this 
sense, the periodic tables shown in Fig. [T] offer 
a simple guide to experimentalists by indicating 
which compositions are expected to be easy and 
which are expected to be hard to synthesize, based 
on thermodynamic considerations. 

Before entering in a detailed discussion on the 
electronic properties of the new stable phases 
found in our simulations, a few general remarks 
are in order, (i) There are some basic atomic ar- 
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Figure 2: Effective hole masses as a function of 
the PBE(+U) band gap for all compounds that are 
within 70meV/atom from the convex hull. The 
color of the dots gives the distance to the hull and 
follows the same scale as in Fig. [TJ 

rangements that are energetically favorable (with 
some exceptions) for many compounds of the 
(Cu,Ag,Au)XC >2 family. The most common ones 
are the delafossite structure of, e.g., CuFeC^ 
(space groups 166 and 194) and the tetragonal 
structure of CuCaC >2 (space group 123), together 
with some slightly distorted, and therefore less 
symmetric, variants. We stress, however, that 
some recurrent crystal structures that we identi¬ 
fied do not have any representative compound in 
databases, (ii) A majority of the structures are lay¬ 
ered, but there is a substantial variety in the ge¬ 
ometry of the layers, (iii) We observe a mixture 
of semiconducting and metallic structures, and for 
some compositions we can even find semiconduct¬ 
ing and metallic phases separated by only a few 
meV/atom. 

Having the crystal structures of new stable and 
metastable (Cu,Ag,Au)X02 compounds opens the 
way for a series of further theoretical studies. In 
fact, we can now calculate a wealth of physi¬ 
cal properties using the numerous methods and 
codes available for theoretical spectroscopy. Our 
main motivation for the choice of this class of sys¬ 
tems is the p-type conductivity measured in Cu 
delafossites. Therefore, we calculated the aver¬ 
age hole effective mass and the electronic band- 
gap for the ground-state of all the stoichiometries 
studied. All calculations were perform with using 
pymatgen® and BoltzTrap 49 software pack- 
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ages. Following the same approach as Ref. 12 we 
calculated the averaged hole effective mass tensor 
for a carrier concentration of 10 18 cm -3 and a tem¬ 
perature of 300 K. We then used the higher limit 
estimation for m* h (see Supplemental Information 
of Ref. fl2l) . The results are summarized in Fig.[2j 
In the figure, the color of the points indicates ther¬ 
modynamical stability, respecting the same color 
scale as in Fig.[TJ We observe that the numbers we 
indicate for the energy gap should be taken with 
care as they are Kohn-Sham gaps obtained with the 
PBE(+U) exchange-correlation functional, and are 
systematically underestimated. 5051 The true gaps 
are substantially larger (usually at least twice as 
large) than the ones we indicate. Moreover, sys¬ 
tems that are metallic within PBE(+U) can some¬ 
times be semiconducting in experiment, however 
the inverse never happens. Note that band disper¬ 
sions, and therefore effective masses, are generally 
less sensitive to the choice of the approximation 
used for the calculation of band structures. 

For p-type transparent conductive materials we 
desire low-hole effective masses (to improve con¬ 
ductivity), and high energy gaps (to ensure trans¬ 
parency). Approximately in the center of Fig. [2] 
we find CuA 1C>2, the compound where p-type elec¬ 
trical conductivity in transparent thin-films was 
discovered for the first time. 31 Its gap is suffi¬ 
ciently large (note that the PBE gap of 1.8 eV is 
substantially lower than the experimental gap of 
around 3.5 eV (see, e.g., Refs. 15015II and refer¬ 
ences therein), but the hole effective mass is rel¬ 
atively high. Indeed, the p-type conductivity in 
CuA 1C>2 is still too low for technological applica¬ 
tions. From Fig. [2] we can see that several materi¬ 
als have the potential to outperform CUAIO 2 , such 
as AgScC> 2 , AuScC> 2 , AuYC> 2 , AgAlC> 2 , AgBiC> 2 , 
etc. Due to the presence of (expensive) noble met¬ 
als (Ag and Au), it is unlikely that such mate¬ 
rials can directly find large scale applications in 
technology, nevertheless an experimental study of 
conductivity in new interesting phases could give 
valuable ideas on if (and eventually how) hole mo¬ 
bilities can be increased beyond present limits. 

In tables [lj [2] and [3] we summarize the infor¬ 
mation concerning the energy distances from the 
convex hull, the band gaps and the hole effective 
masses. 

These results call for more detailed experimen- 



Figure 3: The delafossite structure (space group 
166). 

tal studies of this restricted set of compounds, that 
should be accompanied by more accurate calcula¬ 
tions of the electronic band gaps (i.e. using GW 
approaches beyond standard density functional 
theory 5051 ) and of possible defects/impurities for 
p-type dopability. From a purely theoretical point 
of view, the results of Fig. ^demonstrate that, with 
our approach, we are now able to go all the way 
from a simple stoichiometry to the estimate of rel¬ 
evant material properties for a vast class of mate¬ 
rials. 

We will now analyze more in detail the struc¬ 
tural and electronic properties of the different sub¬ 
classes of compounds that we identified. We will 
extend this discussion to phases that are closer than 
50 meV to the convex hull. 

Delafossites 

The trigonal (space group 166) or hexagonal 
(space group 194) delafossite structures (depend¬ 
ing on the stacking sequence) are the crystal struc¬ 
ture assumed by all stable compounds with X be¬ 
longing to group IIIA (Al, Ga, In), IIIB (Sc, Y, 
La), or X=Cr, Fe, Co, Ni, and Rh. Also AgT102 
and AUTIO 2 (but not CUTIO 2 ) present a low- 
energy delafossite structure. This phase is char¬ 
acterized by XO 2 planes separated by flat hexag¬ 
onal (Cu,Ag,Au) planes. The energy differences 
between the trigonal and hexagonal phases is al¬ 
ways very small (of the order of few meV/atom), 
sometimes smaller than the precision in our cal¬ 
culations. Thermodynamically stable delafossite 
structures not present in databases are AgT 102 , 
AuCr02, AuFe02, AUC 0 O 2 , AuNi02, AUTIO 2 
and AuLa 02 - 
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(a) (b) 



Figure 4: The structures of (a) CuHC >2 (space 
group 31), (b) AgHC >2 (space group 8), and 
(c) AuHC >2 (space group 2). 


(a) (b) 


Figure 5: The structures of (a) CuLi02 (space 
group 2) and (b) (Ag,Au)Li02 (space group 53). 

Group IA 

Hydrogen We found a series of distinct struc¬ 
tures for (Cu,Ag,Au)H02 (see Fig. [4]). The 
Cu compound crystallizes in a layered structure 
(space group 31) composed of zigzag layers of 
CuHC> 2 . It is just above the convex hull at 
8 meV/atom, and is an indirect-gap semiconduc¬ 
tor with a PBE gap of 0.5 eV. The Ag compound 
shows flat layers of Ag (in a square lattice) sepa¬ 
rated by zigzag lines of OH and O 2 dimers. This is 
a monoclinic (space group 8), metallic, thermody¬ 
namically stable phase. Finally, AUHO 2 is a low- 
symmetry structure (space group 2), composed by 
flat Au layers separated by aligned HO 2 units. It 
is 11 meV/atom above the convex hull and it is 
an indirect gap semiconductor with a PBE gap of 
0.7 eV. 

Lithium The structures of (Cu,Ag,Au)Li02 are 
composed of lines of alternating (Cu,Ag,Au) and 



Figure 6: The structure of 

(Cu,Ag,Au)(K,Rb,Cs)02 (space group 63). 

Li bonded by O in a three-dimensional arrange¬ 
ment (see Fig. [5]). The Ag and Au compounds 
crystallize in an orthorhombic structure (space 
group 53). A distortion of this structure leads 
to the monoclinic phase (space group 12) found 
in the databases for CuLiC> 2 , and a large distor¬ 
tion of this structure reduces the symmetry to the 
space group 2 that we find as the ground-state of 
CuLiC> 2 . All compounds are indirect band-gap 
semiconductors, with PBE gaps of 0.5 eV (Cu), 
0.5 eV (Ag), and l.leV(Au). 

Sodium The compounds containing Na result in 
two structures very close in energy and competing 
for the ground-state: the orthorhombic structure 
(see Fig. [6]> that is the structure of the ternary ox¬ 
ides containing K, Rb, and Cs, and the monoclinic 
phase that is a slightly distorted version of the 
structure of (Ag, AujLiCh. For CuNaC >2 we find 
that the former is the ground state (a mere 7 meV 
per atom below the monoclinic structure present 
in databases), while for Ag and Au the latter is 
more stable. These compounds are all indirect- 
gap semiconductors with PBE gaps of 0.4 eV (Cu), 
0.6 eV (Ag), and 1.1 eV (Au). 

Potassium, Rubidium, Cesium The com¬ 
pounds of the form (Cu,Ag,Au)(K,Rb,Cs)02 crys¬ 
tallize in an orthorhombic structure (space group 
63), characterized by corrugated square planes of 
(K,Rb,Cs) intercalated with flat stripes of CUO 2 
where each Cu atom shares four O atoms with two 
neighboring Cu (see Fig. [6]). All these materials 
are semiconducting with an indirect PBE gap be¬ 
tween 0.8 and 1.5 eV that increases with the size 
of the atoms from Cu to Au and from K to Cs. 
Note that this is also the lowest energy structure 
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(a) 


(b) 


Figure 7: The structures of (a) CuMgC >2 (space 
group 67) and (b) AuMgC >2 (space group 65) 

(a) (b) 


(c) 



Figure 8: The structures of (a) CuCa02 (space 
group 123), (b) AgCa02 (space group 12), and 
(c) CuSr02 (space group 63). 

we found for CuNaOo. 

Group IIA 

Magnesium The element Mg is contained in 
two ternary oxides with Cu (43 meV/atom from 
the convex hull) and with Au (36 meV/atom). 
They are both orthorhombic, but quite dissimilar: 
the structure CuMgC^ (space group 67) is a dis¬ 
torted variation of the tetragonal unit-cell (space 
group 123, see Fig. [8]) of, e.g., CuCaC> 2 . AuMgC >2 
is more complicated: its crystal structure (space 
group 65) is similar to the one of (Ag,Au)LiC >2 
(see right panel of Fig. [5]). Both structures are 
metallic. 

Calcium, Strontium, Barium Most of the 
(Cu,Ag,Au)(Ca,Sr,Ba)02 ternaries crystallize in a 
tetragonal lattice with space group 123 (CuCaC> 2 , 
AgBaC> 2 ) or in a distorted version with a mono¬ 
clinic space group 12 (AgCaC> 2 , AgSrC^, AuSrC> 2 , 
AuBaC> 2 ). These two lattices are represented 


in the top panels of Fig. [8j On the other hand, 
CuSr 02 and CuBa 02 have an orthorhombic struc¬ 
ture (space group 63) that resembles to some to 
extent the one of AuMg02. Most of these mate¬ 
rials are metallic, with the exception of the ones 
that crystallize in the monoclinic phase: AgSr 02 
has a PBE gap of 0.6 eV, AuSr02 of 1.4 eV and 
AuBa02 of 1.2 eV. AgCa02 has zero PBE gap, 
but the small density of states at the Fermi surface 
and the results for the other compounds suggest 
that this may just be due to the well known gap 
underestimation of the PBE, and that a quasipar¬ 
ticle gap may open up when more sophisticated 
methods are used. Concerning the thermodynam¬ 
ical stability of these compounds, CuCaC >2 is at 
3 meV/atom above the hull, while AuSrC >2 is at 
33 meV and AuBaC >2 at 31 meV. 

Transition metals 

Vanadium The monoclinic structure (space 
group 12) of CuVC >2 is very similar to the one 
of AgPb02 (see Fig. |T2j). This compound is has 
a gap of 1.0 eV and is 11 meV/atom above the 
convex hull. 

Manganese The oxide phases containing Mn 
(space group 12) can be seen as distorted de- 
lafossites, where the O atoms are slightly dis¬ 
placed from their symmetry positions (CuMnC> 2 , 
AuMnC> 2 ), or with a different stacking (AgMn02). 
The Au compound is metallic whereas the Ag and 
Cu have a PBE gap of 0.4 and 0.1 eV. 

Palladium, Platinum The two Cu compounds 
with Pd and Pt crystallize in the orthorhom¬ 
bic structure of (Ag,Au)LiC >2 (see right panel of 
Fig. [5]). While the Pd compound is thermodynami¬ 
cally stable, CuPt0 2 is 6 meV/atom above the con¬ 
vex hull. Both are metallic. On the other hand, the 
two Ag phases have a monoclinic structure (space 
group 10) similar to AuPb02 (see Fig. |T2]>, but 
where the Ag planes are not distorted (increasing 
the symmetry of the system). AgPd 02 is stable 
with a small PBE gap of 0.1 eV. AgPt02 has an in¬ 
direct gap of 0.3 eV, lying 19meV above the hull. 
Finally, AuPd02 it is a metallic distorted delafos- 
site structure (space group 12) similar to the one 
of (Cu,Ag,Au)Bi02 (see Fig. [13] ). 
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Figure 9: Structure of CuAgC >2 (space group 47). 



Figure 11: The structure of CuT 1C>2 (space group 
11 ). 


Copper, Silver, Gold CuAgC >2 has a orthorhom¬ 
bic lattice characterized by square flat layers of Ag 
intercalated with ribbons of C 11 O 2 . The main dif¬ 
ference to the structure of CuPb02 (see Fig. |T2]) is 
that the four O atoms are shared in this case with 
two Cu neighbors and not with four. On the other 
hand, CuAu 02 crystallizes in a monoclinic struc¬ 
ture (space group 12) very similar to the Au de- 
lafossite structure, with Au hexagonal flat planes 
separated by Cu 02 layers, but where the O atoms 
are slightly displaced thereby reducing the sym¬ 
metry from trigonal to monoclinic. This is, in fact, 
the structure that can be found in the databases for 
CuAg02. Finally, AgAu02 can again be seen as 
a distorted Au delafossite, leading to a monoclinic 
phase with space group 14. While CuAg02 and 
CuAu 02 are metallic, AgAu 02 is an indirect gap 
semiconductor with a PBE gap of 0.6 eV. 



Figure 10: Structure of (a) AgCd02 (space group 
12) and (b) AgHg 02 (space group 2). 


Cadmium, Mercury CuHg02 has the same 
monoclinic structure as CuAu 02 , and so it can be 
seen as a deformed Hg delafossite. It is a metal and 
it lies 43 meV/atom above the convex hull of ther¬ 
modynamic stability. On the other hand, AgCd02 
crystallizes in a deformed version (see Fig. [TO]) 
of the tetragonal structure of, e.g., CuCa 02 (see 
Fig. [8]). This is a metal, 23 meV/atom above the 


hull. The next compound, AgHg02 has as ground 
state a low symmetry structure (space group 2, 
see right panel of Fig. ra that can be seen as 
a small distortion of the orthorhombic phase of 
(Ag,Au)Li02 (see right panel of [5]>. This is a 
metal, 27 meV above the hull. Finally, AuCd02 
has a similar crystal structure as AgHg 02 , also 
with space group 2, but with the metal atoms form¬ 
ing lines of Au and of Cd (i.e., they are not alter¬ 
nating). This is a quasi-direct gap semiconductor 
with a PBE gap of 0.9 eV. 

Thallium The ground-state of CuT 1C>2, in con¬ 
trast to the Ag and Au compounds that crystal¬ 
lize in the delafossite structure, is a monoclinic 
lattice (space group 11), characterized by hexag¬ 
onal channels made of T1 filled with flat stripes of 
Cu0 2 . This phase is an indirect band-gap semi¬ 
conductor with a PBE gap of 0.4 eV. A (metal¬ 
lic) delafossite structure is also present as a meta¬ 
stable phase, around 14 meV/atom higher than the 
ground-state. 

Lead (Cu,Ag,Au)Pb02 oxides are characterized 
by Pb planes separated by (Cu,Ag,Au)C >2 layers 
(see Fig. 0. For CuPb02, Pb forms a square 
lattice while each Cu shares four O atoms with 
four Cu neighbors, while for the Ag compound, 
Pb forms a hexagonal lattice with isolated AgC >2 
units in between. Finally, the AuPb02 structure 
is somewhat intermediate between the Cu and Ag 
compounds: the Pb layers form a distorted hexag¬ 
onal lattice and each Au has two O atoms and 
shares an extra two. The Cu phase (space group 
74) is metallic and 43 meV above the convex hull, 
the Ag compound (space group 12) is a 0.3 eV 
semiconductor and 32 meV above the hull, and 
finally Au forms a semiconducting phase (space 
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(a) (b) 



Figure 12: Structures of (a) CuPbC >2 (space 
group 74), (b) AgPbC >2 (space group 12), and 
(c) AuPbC >2 (space group 2). 



Figure 13: Distorted delafossite structure of 
(Cu,Ag,Au)BiC >2 (compare with Fig. [3j). 

group 2) with an indirect PBE gap of 0.4 eV and it 
is 15 meV above the convex hull. 

Bismuth The ternary oxides with Bi are 
indirect-gap semiconducting phases with PBE 
gaps of 1.0 eV (Cu), 1.4 eV (Ag), and 1.3 eV 
(Au). Their monoclinic space group 11 corre¬ 
sponds to an atomic arrangement that can be seen 
as a distorted delafossite, with a dimerization of 
the chains of Ag and BiC >2 (see Fig. [13]). The de¬ 
lafossite structure remains, however, a meta-stable 
phase of this composition, around 40meV/atom 
higher than the ground-state. Although not strictly 
stable thermodynamically, these structures are re¬ 
markably close to the convex hull at 8 meV/atom 
(Cu), 3 meV/atom (Ag), and 40 meV/atom (Au). 

Halogens 

Fluorine (Cu,Ag,Au)FC >2 compounds crystal¬ 
lize in a monoclinic lattice (space group 4) char¬ 



Figure 14: The structure of (Cu,Ag,Au)F02 (space 
group 4). 


(a) (b) 



Figure 15: The structure of (a) Ag(Cl,Br)02 
(space group 1), (b) Au(Cl,Br)C >2 (space group 
15), and (c) Cu(Cl,Br)02 (space group 1). 


acterized by (Cu,Ag,Au)F02 layers (see Fig. |T4j). 
These are all semiconducting structures, with 
PBE electronic gaps of 0.8 eV (Q 1 FO 2 ), 0.6 eV 
(AgF0 2 ), and 1.1 eV (AuF0 2 ). 

Chlorine, Bromine Cu(Cl,Br)C >2 compounds 
crystallize in a low symmetry triclinic (space 
group 1) lattice (see panel c of Fig. |T5]). This 
has some similarities to the monoclinic structure 
of (Cu,Ag,Au)FC >2 (see Fig. [14]), in the sense that 
they share the same CUO 2 subunits. However, 
in the latter structure each Cu shares four F with 
four neighboring Cu atoms, while in the former 
each Cu only shares two (Cl,Br) with two other 
Cu. These are thermodynamically stable struc¬ 
tures, that are indirect gap semiconductors with 
PBE gap of 0.8 (CUCIO 2 ) and 0.9 (CuBrC> 2 ). 

Ag(Cl,Br)02 structures are fundamentally dif¬ 
ferent from all other we found in our study. It is 
composed of flat hexagonal layers of Ag(Cl,Br), 
separated by a layer of O 2 molecules (see panel a 
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Figure 18: The structure of AgAsC >2 (space group 
5). 



Figure 16: The structure of AgBC >2 (space group 
9). 


(a) (b) 



Figure 17: Two phases of AgNC> 2 , (a) the 
ground-state monoclinic (space group 5)and an or¬ 
thorhombic (space group 44) structures. 


group 5), lying 37 meV/atom above the hull. It 
is composed of lines of alternating Ag and NO 2 . 
This is a deformation of the orthorhombic struc¬ 
ture (space group 44, see Fig. E3 that can be 
found in the databases, and that we found to 
be 5 meV/atom above the ground state. Both 
structures are indirect gap semiconductors, with 
PBE gaps of 1.4 eV (monoclinic) and 1.8 eV (or¬ 
thorhombic). 


of Fig. [T5| >, forming a low-symmetry triclinic lat¬ 
tice (space group 1). These turn out to be semi¬ 
conductors with PBE gaps of 0.6 eV (AgCICM and 
0.7 eV (AgBr0 2 ). 

Au(Cl,Br)02 compounds crystallize in a lattice 
that is composed of Au(Cl,Br) layers intercalated 
with O 2 molecules (see panel b of Fig. [15]) . How¬ 
ever, in this case the Au(Cl,Br) are not flat, but 
form zigzag stripes with each Au bonded to two 
(Cl,Br) and vice-versa, leading to a monoclinic lat¬ 
tice (space group 15). AuC 1C>2 is a semiconductor 
with an indirect gap of 0.3 eV, while AuBrC >2 turns 
out to be a metal with the PBE. 


Others 

Boron AgB02 has a monoclinic structure (space 
group 9) composed by slightly buckled Ag layers, 
separated by BO 2 chains, where each B is con¬ 
nected to three O, two of which are shared with 
two adjacent B atoms (see Fig. |T6j). Note that this 
is the same structure we had previously found for 
CUBO 2. 28 AgB02 is 33 meV/atom above the hull, 
and it is a semiconductor with a PBE quasi-direct 
gap of 1.4 eV. 

Nitrogen The lowest energy phase of AgNC >2 
that we found is a monoclinic structure (space 


Arsenic Also for AgAs02 we found a mono¬ 
clinic structure (space group 5, see Fig. p~8]), com¬ 
posed of corrugated planes of AgAsOo. This phase 
is 38 meV/atom above the hull, with a PBE indi¬ 
rect gap of 2.4 eV. 


Conclusion 

In summary, we performed structural prediction 
runs for 183 oxide phases containing Cu, Ag, and 
Au. From our runs we predict that there are 81 
thermodynamically stable or quasi-stable (within 
20meV) compositions, out of which only 36 are 
included in available databases. These numbers 
should be compared to the 3 systems that were 
found in a previous study using high-throughput 
techniques combined with machine learning. Nev¬ 
ertheless, we believe that structural prediction 
should not be seen as a competitor of more tra¬ 
ditional high-throughput techniques. In fact, both 
methods complement each other, and both should 
be used in synergy to speed up the experimental 
process of materials discovery. We tried to do the 
first steps in this direction, by combining structural 
prediction and a search based on prototype struc¬ 
tures. 

The subset of oxides that we explored include 
delafossite CUAIO 2 , the first delafossite p-type 
transparent conductive oxide. By calculating the 
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band gaps and hole effective masses of the new sta¬ 
ble compounds that we identified, we could reveal 
some potential candidates to outperform CuAlC^ 
as p-type transparent conductor. These few com¬ 
pounds deserve in our mind further consideration, 
both from an experimental point of view (e.g., syn¬ 
thesis, structural and electronic characterization) 
and a theoretical point of view (e.g., study of dopa- 
bility, more accurate band structure and transport 
calculations). 

Certainly, structural prediction is complex and 
numerically expensive. However, with our results 
we aim at demonstrating that, using modern super¬ 
computers, it is possible to use such techniques to 
investigate a large number of chemical composi¬ 
tions, paving the way for a large scale search for 
new materials. 
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Table 1: Space groups (Spg), distance to the con¬ 
vex hull (E h u n), gaps (Gap) and average hole ef¬ 
fective mass (mf) for semiconducting structures of 
the form CUXO 2 lying below 50 meV/atom from 
the convex hull. 


Structure 

E hull 

Spg 

Gap 

m h 

CuH0 2 

8 

31 

0.5 

2.23 

CuFi02 

0 

2 

0.5 

17.9 

CuF0 2 

0 

4 

0.8 

11.8 

CuNa0 2 

0 

63 

0.4 

3.04 

CuMgOo 

43 

67 

0.0 

- 

CuA10 2 

0 

166 

1.8 

2.66 

CuC10 2 

0 

1 

0.8 

48.9 

CuK0 2 

0 

63 

0.8 

13.7 

CuCa02 

3 

123 

0.0 

- 

CuSc0 2 

0 

194 

2.4 

4.50 

CuV0 2 

11 

12 

1.0 

3.42 

CuCr0 2 

0 

166 

1.6 

5.52 

CuMn0 2 

0 

12 

0.1 

0.06 

CuFe02 

1 

166 

0.9 

2.02 

CuCo0 2 

0 

166 

1.1 

5.63 

CuGa0 2 

0 

166 

0.8 

2.23 

CuBr0 2 

0 

1 

0.9 

9.67 

CuRb0 2 

0 

63 

0.8 

8.98 

CuSr02 

0 

63 

0.0 


CuY0 2 

0 

194 

2.6 

4.01 

CuRh0 2 

0 

166 

0.7 

2.75 

CuPd0 2 

0 

53 

0.0 

- 

CuAg0 2 

0 

47 

0.0 

- 

Culn0 2 

7 

166 

0.3 

0.57 

CuCs02 

0 

63 

0.9 

5.24 

CuBaOo 

0 

63 

0.0 


CuLa02 

2 

166 

2.7 

4.04 

CuYb0 2 

0 

123 

0.0 

- 

CuPt0 2 

6 

53 

0.0 

- 

CuAu0 2 

0 

12 

0.0 


CuHg0 2 

43 

12 

0.0 


CuT10 2 

31 

11 

0.4 

4.40 

CuPb0 2 

43 

74 

0.0 


CuBi02 

8 

11 

1.0 

4.40 
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Table 2: Space groups (Spg), distance to the con¬ 
vex hull (E h u n), gaps (Gap) and average hole ef¬ 
fective mass (mp for semiconducting structures of 
the form AgXC >2 lying below 50 meV/atom from 
the convex hull. 


Structure 

E hull 

Spg 

Gap 

m h 

AgH0 2 

0 

8 

0.0 

— 

AgLi0 2 

0 

53 

0.5 

1.85 

AgB02 

33 

9 

1.4 

1.61 

AgN0 2 

37 

5 

1.4 

2.57 

AgF02 

0 

4 

0.6 

7.61 

AgNa0 2 

0 

12 

0.6 

4.79 

AgA10 2 

0 

166 

1.4 

1.49 

AgC10 2 

0 

1 

0.6 

2902 

AgK0 2 

0 

63 

0.9 

5.80 

AgCa0 2 

0 

12 

0.0 

- 

AgSc0 2 

0 

194 

2.1 

2.58 

AgCr0 2 

0 

166 

1.7 

2.48 

AgMn02 

0 

12 

0.4 

1.33 

AgFe02 

0 

166 

1.1 

0.90 

AgCo02 

0 

194 

1.2 

3.68 

AgNi0 2 

20 

166 

0.0 

- 

AgGa0 2 

0 

166 

0.6 

1.09 

AgAs0 2 

38 

5 

2.4 

3.25 

AgBr02 

0 

1 

0.7 

450 

AgRb0 2 

0 

63 

1.0 

8.80 

AgSr0 2 

0 

12 

0.6 

1.27 

AgY0 2 

0 

194 

2.4 

3.82 

AgRh0 2 

38 

166 

0.5 

2.14 

AgPd0 2 

0 

10 

0.1 

0.03 

AgCd0 2 

23 

12 

0.0 

- 

Agln02 

0 

166 

0.2 

0.43 

AgCs0 2 

0 

63 

1.1 

4.54 

AgBa0 2 

0 

123 

0.0 

- 

AgLa02 

0 

194 

2.8 

3.85 

AgPt0 2 

19 

10 

0.3 

3.30 

AgAu02 

0 

14 

0.6 

1.80 

AgHg0 2 

27 

2 

0.0 

- 

AgT10 2 

0 

166 

0.0 

- 

AgPb0 2 

34 

12 

0.3 

0.03 

AgBi02 

3 

11 

1.4 

2.38 


Table 3: Space groups (Spg), distance to the con¬ 
vex hull (E hull), gaps (Gap) and average hole ef¬ 
fective mass (mp for semiconducting structures of 
the form AuX02 lying below 50 meV/atom from 
the convex hull. 


Structure 

E hull 

spg 

Gap 

m h 

AuH0 2 

12 

2 

0.7 

1.68 

AuLi0 2 

0 

53 

1.1 

3.78 

AuF0 2 

0 

4 

1.1 

3.97 

AuNa0 2 

0 

12 

1.1 

7.72 

AuMg0 2 

37 

65 

0.0 


AuA10 2 

12 

166 

0.6 

0.71 

AuC10 2 

0 

15 

0.3 

0.96 

AuK0 2 

0 

63 

1.4 

9.20 

AuSc0 2 

14 

166 

1.8 

0.63 

AuCr0 2 

0 

166 

1.2 

1.65 

AuMn0 2 

0 

12 

0.0 

- 

AuFe0 2 

0 

166 

0.6 

0.59 

AuCo0 2 

0 

194 

0.5 

1.03 

AuNi0 2 

0 

194 

0.0 

- 

AuGa0 2 

6 

166 

0.0 

- 

AuBr0 2 

0 

15 

0.0 

- 

AuRb0 2 

0 

63 

1.4 

8.45 

AuSr0 2 

33 

12 

1.3 

1.80 

AuY0 2 

3 

194 

2.7 

2.49 

AuPd0 2 

2 

12 

0.0 


AuCd0 2 

0 

2 

0.9 

1.70 

Auln0 2 

18 

166 

0.0 

- 

AuCs0 2 

0 

63 

1.5 

6.59 

AuBa0 2 

31 

12 

1.2 

1.44 

AuLa0 2 

21 

194 

2.7 

4.45 

AuT10 2 

13 

166 

0.0 


AuPb0 2 

15 

2 

0.4 

1. It 

AuBi0 2 

40 

11 

1.3 

1.54 
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The delafossite crystal structure. 
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